Submitted to MNRAS on 04/09/2012 



On the Nature of the TeV Emission from the Supernova Remnant SN 1006 

Miguel Araya 1 & Francisco Frutos 

Space Research Centre ( CINESPA ) 
University of Costa Rica 
San Jose 2060, Costa Rica 

ABSTRACT 

We present a model for the non-thermal emission from the historical supernova 
remnant SN 1006. We constrain the synchrotron parameters of the model with archival 
radio and hard X-ray data. Our stationary emission model includes two populations 
of electrons, which is justified by multi-frequency images of the object. From the set 
of parameters that predict the correct synchrotron flux we select those which are able 
to account, either partly or entirely, for the gamma-ray emission of the source as seen 
by HESS. We use the results from this model as well as the latest constraints imposed 
by the Fermi observatory and conclude that the TeV emission cannot be accounted for 
by 7T° decay from high-energy ions with a single power-law distribution, of the form 
dN pmt JdE p oc E p s , and s > 2. 

Subject headings: gamma-ray: observations; ISM: supernova remnants; 
ISM:individuals:SN 1006, acceleration of particles, radiation mechanisms: non-thermal 

1. Introduction 

Supernova remnants (SNRs) are thought to be the main source of galactic cosmic rays, high- 
energy (E p < 10 15 eV) particles (mostly protons) that populate the galaxy, and whose distri- 
bution is a power-law. SNRs are energetically favorable and the Fermi mechanism is the natu- 
ral process by which these objects are able to inject part of their energy into charged particles 



(Bel l 19781 : Blandford & Eichler 1987f) . Such diffusive shock acceleration (DSA) could account 



for the observed shape of the cosmic ray distribution. 
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A direct detection of high-energy protons and nuclei at the sources is hard to obtain, however, 
and the preferable mechanism to accomplish this task is the detection of MeV-TeV 7-ray emission 
from the decay of neutral pions resulting from inelastic collisions between high-energy ions and 
ambient ions (the so-called hadronic 7-ray emission). This hadronic signal then has to be separated 
from (mostly MeV) radiation from the galactic background as well as from (MeV-TeV) radiation 
from high-energy electrons in SNRs (the so-called leptonic 7-ray emission). 

High-energy electrons were first detected in SNRs through their non-thermal radio spectrum. 
This radiation is produced when accelerated electrons move in the presence of magnetic fields, 
which result either from shock compression of the interstellar field or from amplification by cos- 
mic ray instabilities (see lSchure et al. 20121 for a recent review). These electrons typically have 
energies in the GeV range. 

The presence of TeV electrons in young SNRs has been confirmed by X-ray observations. 
High-resolution observations of young SNRs (e.g., taken by the Chandra X-Ray Observatory) 
show synchrotron emission associated with the forward shock (e.g., 



Gotthelf etal. 2001 



Berezhko et al. 2002 ; 



Hwang et al. 2002 Rho et al. 2002l : lLong et al. 2003l : IVink & Laming 20031: Berezhko & Volk 2004 : 
Arava etal. 2010). 



The 7-ray leptonic emission results from inverse-Compton (IC) up-scattering of background 
photon fields by the synchrotron-emitting electrons and from bremsstrahlung emission from colli- 
sions between these electrons and ambient plasma. Leptonic and hadronic high-energy emission 
could in principle be distinguished based on their relative flux normalization as well as the intrinsic 
differences in their spectra. However, due to the limited knowledge on the physical parameters of 
some SNRs, the model degeneracy and limited observations, the identification of hadronic signa- 
tures in their non-thermal spectra can be quite challenging. 

Spectral energy distributions (SEDs) of SNRs interacting with molecular clouds (which are 
typically old, with an age > 10 3 years), such as W28, W41, W49B, W51C, IC 443, G8.7-0.1 and 



G359. 1-0.5, seem t o favor a hadronic origin for the emission (lAbdo et al. 2010bl:lM ehault et al. 2011 
Abdo et al. 20~10cl : lAbdo et al. 20091 lAbdo et al. 2010al : kjello et al. 2012l : lHui et al. 201 lb . while" 



other recent studies show that the SEDs of some you ng SNRs with a h ard GeV spectrum are 
probably leptonic-d ominated, such as RX J17 13.7-3946 (lAbdo et al. 2011 ah and RX J0852.0-4622 
(|Tanaka et al. 201 lb . and have spectral slopes at GeV energies that are consistent with the corre- 
sponding IC spectra expected from synchrotron-emitting electrons. 



Other young objects, such as the f amous historical S NRs Cassiopeia A (lAraya & Cui 2010 : 
Abdo et al. 2010db and Tycho's SNR d Arava et al. 201 1|) show a soft GeV spectrum which, it 
might be argued, favors the hadronic scenario (see lAtoyan & Dermer 201 1[ however, for an al- 
ternative explanation for the emission from Tycho's SNR). Recently, an attempt to unify the prop- 
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erties of 7-ray emitting SNRs showed that as the den sity of the envir onment of the SNR increases, 



so does the brightness of the hadronic component (|Yuan et al. 2012|) . 



Regardless of the nature of the high-energy photon emission, it seems clear that the spectral 
shape below and around hundreds of MeVs is fundamental for its identification. The recently- 



launched Fermi satellite (|Atwood et al. 2 009) has allowed probing this region of the spectrum, 
although the increasing PSF of its Large Area Telescope (LAT) at these energies and the high 
galactic background levels may complicate the analysis for some sources. 

Here we present a model for the non-thermal emission from SN 1006, the remnant of a type 



la supernova reported by Chinese and Arab astronomers in 1006 (Ste phenson & Green 2002|) . The 
SNR is located far from the galactic center (galactic lati tude ~ 14.5 °), in a "clean" environment 
with l ow density (~ 0.085 cm" 3 , see liCatsuda et al. 2009 ). at an estimated distance of 2.2 kpc from 



Earth ([Winkler et al. 2003|) . We adopt these values here and develop a two-zone leptonic model 
to find the sets of reasonable parameters that account for the synchrotron emission observed and 
are also consistent with the 7-ray flux as seen by the H.E.S.S. atmospheric Cherenkov telescopes 



( Aceroetal. 2010 ). 



As fundamental part of the data used, we included the latest constraints on GeV observations 
of the source, obtained through analysis of 3.5 years of Fermi LAT data. Although the analysis 
shows that the source is not yet significantly detected at GeV energies, the established upper limit 
for the emission is used to constrain the SED models. 

The results of our simplified model indicate that the TeV flux can only be accounted for 
by leptonic emission from high-energy electrons in a relatively low (~ 30 /xG) magnetic field, 
hadronic emission from a hard power-law particle distribution, or by a mixed leptonic-hadronic 
scenario. In the later cases, hadronic emission from ions with a soft power-law distribution with 
index higher than 2 can be ruled out as the main cause of TeV photons. 



2. Data 

The SED of SN 1006 contains archival radio observations from different instruments, and 



were taken from Dyer et al. (2009). The radio flux density can be described by a power-law with 
an index of ~ —0.57, which can be accounted for by optically-thin synchrotron radiation from a 
particle distribution with index —2.14. 

The X-ray fluxes of the source were obtained from observations taken by the Proportional 
Counter Array (PC A) onboard RXTE on 1996 February. The procedure followed for reducing the 
data is standard for PCA data, the HEASOFT package version 6. 1 1. 1 was used. Only events from 



-4- 



the top layer of PCU2 were used for analysis in the interval from 3.5 keV to 25 keV, and good 
time intervals were generated by filtering out events during South Atlantic Anomaly passages, 
times when SN 1006 was less than 10° above the limb of the Earth and the pointing direction of 
the detectors is more than 0.°02 from the nominal pointing direction. The PCA background spec- 
trum for SN 1006 was calculated using version 3.8 of the FTOOL pcabackestQ and the back- 
ground model pcaJbkgd-cmfaintl7-eMv20051128.mdl. The background- subtracted spectra from 
the source for each independent observation were combined with the FTOOL addspec and ap- 
propriate responses were generated with addrmf . 

The resulting X-ray spectrum was binned to contain 1000 counts per bin and fitted in XSPEC 
version 12.7.c| with an absorbed power-law model after adding 1% systematic error. The column 
density cannot be co nstrained and is fixed at N H = 7 x 10 20 cm~ 2 , according to the literature 



(|Kalberla et al. 20051 . and also see the FTOOL nh which gives archival galactic column density 



values). 

The fit obtained has a reduced x 2 of 1 .08, and the resulting values for the photon index and un- 
absorbed flux (with 90%-confidence errors) are 2.83±g;];| and (2.08 ± 0.11) x 10~ n ergs/cm 2 /s (in 
the energy band 3.5 — 25 keV), respectively. The resulting unfolded sp ectrum, corrected for absorp- 
tion using standard photoelectric cross section values below 10 keV ( Morrison & McCammon 198 3 1, 
is shown in Fig. [Qwith the rest of the SED. 



TeV 7-ray observations of SN 1006 published by the H.E.S.S. Collaboration (|Acero et al. 2010)) 
and the latest constraints imposed by the Fermi LAT obtained here complete the high-energy SED 
of the source. 



2.1. Analysis of Fermi LAT data 

We analyzed the Fermi LAT data gathered between 04 August 2008 and 16 January 2012. 
The analysis was performed with the LAT Science Tools package □ version v9r23pl. The cuts and 
selection criteria applied are those recommended by the Fermi Science Support Center for binned 
likelihood analysis^, which is generally used for large data sets such as this one. We selected diffuse 



For more information about FTOOLS, see http://heasarc.gsfc.nasa.gov/docs/software/ftools . 



For more information about XSPEC, see http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 



3 The spectral absorption factor is e n h°(e)^ where a{E) is the cross section for a photon energy E. 
4 See http://fermi.gsfc.nasa.gov/ssc 



5 See http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/binnedJikelihood_tutorial.html . 



-5- 



events (see lAtwood et al. 2009|) between 500 MeV and 100 GeV within a region of interest (Rol) 
of 30° around the position of the source, RA (J2000)= 15 h 02 m 22. I s , Dec (J2000)= -42°5'49 // , and 
a maximum zenith angle cut of 100°, in order to avoid the background from the limb of the Earth, 
with the tool gt select. Appropriate livetime calculation with the tool gtmkt ime is performed, 
excluding time intervals where the zenith cut intersects the Rol. 

In order to account for all the photons from background sources, an 'xml' file is created which 
conta ins information abo ut source positions an spectra as published in the LAT 2-year Source Cat- 



alog (lAbdo et al. 201 lb). The resulting file contains 260 point sources and 2 extended sources 
(known as Cen A and MSH 15-52), which are modeled according to the emission templates pro- 
vided by the Fermi Collaboration^. Modeling such a large region of the sky is necessary due to the 
large PSF of the telescope. The size of the region and the time span of the observation make the 
analysis computationally extensive. However, the source is located relatively far from the galactic 
center which reduces the uncertainties in the galactic background estimation. 

The current instrument response functions (P7SOURCEJV6) are used throughout the anal- 
ysis, as well as the latest galactic and extragalactic diffuse background components (as specified 
in the files galJ2yearp7v6_v0.fits and isojp7v6source.txt, respectively). The spectral parameters of 
the catalogued sources beyond 30 degrees from the position of SN 1006 are kept frozen during the 
fit, which is performed several times with the use of the optimizer NEWMINUIT until convergence 
is achieved. 



The diameter of the X-ray shell of the remnant, as seen for instance by Chandra (|Long et al. 20031) . 



is about 30 arcmin. The averaged 9 5%-containment ang le for (front+back) events for the LAT is 



about 0.8° at the highest energies (|Atwood et al. 20091) where the resolution is also the highest. 



This would make SN 1006 hard to resolve by the Fermi LAI0. A point source hypothesis for SN 
1006 was then used to derive the main flux constraints presented in this paper. 



included within the software distribution, or downloaded at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yrjcatalog/gll_psc_v07_templates.tgz . 

7 An extended source hypothesis was tested in the likelihood anal ysis as well, follo wing the emission profile for the 
source reported from the Sydney University Molonglo Sky Survey dBock et al. 19991) at 843 MHz, and the results do 
not change significantly with respect to the point-source model. 
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2.1.1. Results of the LAT data analysis 

The emission from SN 1006 at GeV energies is not significant. The resulting value of the 
test-statistics§ (TS) for a point-source at the position of SN 1006 with free spectral parameters 
is only 7.5, which corresponds to a detection significance of only about \/TS ~ 2.7a above the 
background. 

If the 7-ray emission from the source is assumed to follow a power-law distribution, diV 7 / <iE 7 oc 
E~ r (in MeV -1 cm~ 2 s _1 , with T the photon index), with T = 2.14, as would be expected from 
emission of hadronic origin from a high-energy proton population with a distribution similar to the 
one observed for the high-energy electrons, dN pioton /dE p oc E p 2 - 14 (with E p the proton energy), 
the TS value decreases to about 4.9, and the estimated 99%-confidence level upper limit on the 
integrated source flux in the energy interval of the analysis (0.5 — 100 GeV) under this assumption 
is 1.32 x 10~ 9 photons/cm 2 /s. 

The TS values obtained for different assumptions on the index of the power-law photon distri- 
bution by a point source at the position of SN 1006 are found to decrease with higher photon index 
(i.e., for softer 7-ray spectra). For example, the TS value obtained for a photon index Y = 2.0 
is 5.7, and for T = 1.6, 7.2. This latter value corresponds to the slope of the photon distribution 
(i.e., in photons MeV -1 cm~ 2 s _1 ) expected for GeV emission resulting from IC scattering by an 
electron population with index 2.14, as the one observed in SN 1006. Table \T\ shows a summary 
of these results. The upper limit obtained for T = 1.6 is shown in the SED (Fig. [D and is used to 
constrain our leptonic model, described in the next Section. 



3. Model 

3.1. Scales of emission regions and particle distributions 

The observed radio and X-ray synchrotron emission is accounted for with a two-population 
stationary model. The diffuse radio electrons, which we denote as Zone 1, are uniformly distributed 
over a volume V\ and higher energy electrons in Zone 2 associated with the X-ray rim near the 
forward shock region, with a volume V 2 . 



From the radio and X-ray images of the source shown by Cassam-Chena'i et al. (|2008i) . we 
estimated that V 2 is about 11% of the volume of the SNR (or V 2 ~ 1.4 x 10 55 cm 3 , for an spherical 



The modeling of Fermi LAT data i nvolves a likelihood fit where the significance of a source is quantified with the 
value of its TS (see Mattox et al. 19 96, for details). 
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SNR of radius 9.5 pc), and the volume fraction of Zone 1 is about 62% of the volume of SN 1006 

(Vi ~ 7.6 x 10 55 cm 3 ). 

Next, the non-thermal emission from these two populations is calculated, using power-law 
particle distributions in the two Zones, as usually done to account for optically thin synchrotron 
radiation. In fact, in order to account for the X-ray data, a broken power-law particle distribution 
is used in Zone 2, and similarly in Zone 1 as required for the high-energy 7-ray data (see below). 

The synchrotron emission from Zone 1 dominates at radio frequencies, while Zone 2 is bright- 
est at X-ray energies. Below the break energy, the power-law index for both particle distributions 
is fixed at 2.14, which is required by the radio spectrum. 

We calculated the resulting 7-ray emission from these two populations. It included inverse 
Compton scattering of cosmic microwave background photons (IC/CMB) and non-thermal bremsstrahlung 
(NB). For the latter we considered interactions of electrons with other electrons, protons, and fully 
ionized helium. We used an abundance for helium of 0. 1 that of hydrogen. The contributions from 
synchrotron self-Compton scattering are entirely negligible. 

The NB flux is proportional to the density of the target plasma, which for Zone 2 can be taken 
as four time s the local ISM density, nf^ 2 = 0.3 cm~ 3 , as is standard for a strong shock (e.g. 



Drury 1983|) . and for Zone 1 it can be estimated as the mass of the ejecta (1.4 M Q ) plus a com- 
parable mass of shocked ISM material distributed in the volume of the SNR, n^ ncl = 0.04 cm -3 . 
These values result in negligible NB emission for all reasonable values of the model parameters. 

The contribution to the 7-ray emission from Zone 2 is negligible for all the reasonable values 
of the magnetic field in this Zone (see Section [3721) . The particle indices beyond the break are 
2.8 for Zone 1 and 2.9 for Zone 2, necessary to reproduce the 7-ray and X-ray observations, 
respectively. It is beyond the scope of this paper to speculate on the nature or cause of the particle 
breaks, and we limit ourselves to point out the need of such spectral distributions if the emission is 
to be accounted for by leptonic processes alone. 

The two particle distributions in our main leptonic emission model (Fig. [B can be reasonably 
well described by the following functions 

2 14 / 2 \ 0.33 

(,X Z 1 (2-06 x 10- GeV cm-) (J—) "" 1 + ( -J—) and 



dEJV v '\lGeVj \ \^ l m e c 




-0.38 



dN e Zone2 , in , „ / E e V 2 ' 14 ( ( E e 

= 1.84 x 10- 10 GeV" 1 cm" 3 — ^— 1 + 6 

dE e dV v 7 V 1GeV / V V7br n m eC 2 , 

where m e is the mass of the electron, and 7^ nel = 7 x 10 5 and 7fc°" c2 = 3 x 10 6 are the break 
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Lorentz factors for Zone 1 and Zone 2 respectively. The corresponding maximum Lorentz factors 
are 5 x 10 7 and 6 x 10 7 . 



3.2. Magnetic field and energy in electrons 



We adopt a value for the magnetic field in the X-ray rim of SN 1006 (Zone 2) co nsistent 
with constraints on the diffusion and acceleration parameters made by Parizot et al. (|2006|) . Their 
calculations, obtained assuming uniform, isotropic turbulence, are based on observed SNR X-ray 
properties such as the synchrotron cutoff energy, the velocity of the expanding shock and the rim 
thickness. 



The values obtained by Parizot et al. (120061) for the field in the downstream region of SN 1006 
are ~ 60 — 100 fiG, depending on the value of the compression ratio and the energy dependence of 
the diffusion coefficient used. They u sed a value for the p rojected width of the rim that is consistent 
with other estimates (I = 0.2 pc, see lBamba et al. 2003|) and a distance to the remnant of 2.2 kpc. 
Although they used a value for the shock speed of 2900 km/s that is a factor of 1.7 lower than 
more recent estimates (fKatsuda et al. 20091) . which according to their formula for the downstream 
magnetic field translates into a field value around 1.2 times larger. 

The value adopted here is 80 /iG, which results in an energy content of W^ 2 = 7.7 x 10 43 
erg for the corresponding electrons. This value for the magnetic fiel d is also consistent with other 
values obtained from models applied to high resolution X-ray data (pvlorlino et al. 20101). althoug h 
considerably lower compared to some nonlinear kinetic models (~ 150 /iG. see lvolk et al. 2005 ). 
We explored the effects of using such a range of values for the field in Zone 2 in our model (see 
below). 

With the parameter values described so far, the resulting IC emission from Zone 2 is small, 
and considerably below the observed H.E.S.S. flux at TeV energies. At this point we consider 
a leptonic scenario for the very high energy emission from SN 1006. If the 7-ray emission is 
attributed to Zone 1, the magnetic field required there should be relatively low, around 10 /iG, 
resulting in an energy content of W^""" 1 = 4.2 x 10 45 erg. 

In such a leptonic scenario, the break energy in the particle distribution in Zone 1 should 
be lower than the corresponding value for Zone 2, in order to account for the observed H.E.S.S. 
emission at the highest energies without increasing the GeV fluxes, which are constraint by the 
Fermi LAT observations. 



The resulting modeled emission as well as the observations are shown in Fig. \T\ As can be 
seen, the exact slope of the H.E.S.S. data is hard to reproduce with the model. This is the case for 
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a wide range of model parameters. 



3.3. On the parameter degeneracy 

If the model parameters for Zone 1 are now fixed to the values described above, the required 
magnetic field in Zone 2 cannot be lower than about 30 /xG, for a corresponding energy in electrons 
in this Zone of 5.6 x 10 44 erg and a higher maximum Lorentz factor of 1.2 x 10 8 . The predicted 
7-ray SED flux from Zone 2 for this field would be about 5 times lower than the emission from 
Zone 1 at most frequencies, but the lower limit on the magnetic field results from fitting the 7-ray 
data with the contributions from both Zones. In other words, this range of values is allowed by the 
uncertainty in the data. 

Conversely, if a higher value of the magnetic field in Zone 1 (the radio emitting electrons) is 
used, the TeV fluxes could instead be attributed to Zone 2, where the field should be lower. It was 
found that for a magnetic field in Zone 1 of 40 /iG, the required field in Zone 2 should be about 
17 /xG. Also, the break and maximum Lorentz factors would be higher than the values used above, 
and both electron populations would have comparable total energies. However, this field in Zone 2 
is not realistic for the observed widths of the filaments. For this reason, we adopted the parameters 
described in Sections 13.1 l and |X2l for the rest of the analysis. We refer to this particular parameter 
selection as our main leptonic model. 



3.4. Hadronic emission 



We studied the effect of adding a component of 7-ray emission from hadronic interactions 
to the modeled non-thermal SED of SN 1006. The hadronic emission is calculated as in Kamae 



et al. (|2006() for power-law particle distributions with various slopes. Figures |2] and |3] show the 
corresponding models obtained for different hadronic emission levels and particle indices. 

When assessing the possible amount of hadronic flux, the parameters for Zone 2 were fixed 
for all the models to the values shown in Sections [3TT1 and [3T2l but the amount of IC emission from 
Zone 1 is varied to maintain consistency with the data as required. 

Figure |2] shows the result of emission models with a steep hadronic population described by 
a single power-law with index 2.14. This value also corresponds to the index of the distributions 
of the leptonic components. Such a value for the index of the proton distribution is typical of 
hadronic 7-ray em ission used in current SED models, bo t h of SNRs and SNRs interacting with 
molecular clouds (krava & Cui 20ip[ kbdo et al. 2010al kbdo et al. 201 Pel : kbdo et al. 201 Pel : 
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Arava et al. 201lUGiordano et al. 20121) . 



For the case of the hadronic-dominated scenario, the value of the magnetic field in Zone 1 is 
40 /iG, resulting in an energy content in electrons of 4.6 x 10 44 erg. The shape of the spectrum 
of the radio-emitting electrons is the same as before, but the required maximum Lorentz factor in 
Zone 1 is lower than in the leptonic-dominated scenario. The maximum particle energy reached 
by protons, as required by the H.E.S.S. data, is 90 TeV. 

As can be seen in Fig. [21 the hadronic model over-predicts the flux at GeV energies and 
therefore is not consistent with the limit imposed by observations. In fact, the same conclusion 
can be obtained for any power-law proton distributions with an index higher than 2.0, due to 
the shape of the resulting spectra (the 7-ray flux from hadronic interactions maps the particle 



spectrum directly, see, e.g., Kamae et al. 2006). Therefore, hadronic TeV emission from a cosmic 
ray population with index greater than 2.0 at all energies at the level observed by H.E.S.S. is 
discarded for SN 1006. 

For a proton distribution with index 2.14, Fig. [2] also shows a scenario where the hadronic 
emission is the maximum possible that is also compatible with observations. This is achieved with 
a magnetic field and energy in Zone 1 of 12 fiG and 3.0 x 10 45 erg, respectively. 

Similar models can be obtained for a different index for the proton distribution. For a proton 
distribution with index 2.0 (the value obtained from DSA calculations for strong shocks in the 
test-particle limit), models are shown in Fig. [3] for the case where the emission is dominated by 
hadronic processes (and which is obtained by suppressing the IC emission with a magnetic field in 
Zone 1 of 50 fiG and an energy in electrons of 3.0 x 10 44 erg) and for a case where the leptonic 
and hadronic fluxes are comparable at TeV energies (with _£? Zonel = 16 /xG and an energy in Zone 
1 of 2.1 x 10 45 erg), which we argue is more consistent with the data (particularly with the LAT 
upper limit). In these two cases, the proton cosmic ray spectrum extends up to 45 TeV. 



4. Discussion and Conclusions 

From the simple model presented here, we can conclude that if a value for the magnetic field 
in Zone 2 is used which would reproduces the widths of X-ray filaments near the forward shock 
of SN 1006 (~ 80 /iG), then a leptonic scenario for the 7-ray emission implies that high-energy 
photons should originate mainly from radio-emitting electrons in Zone 1. In this scenario, the 
required electron spectrum is described by a broken power-law with indices 2.14 (as required by 
the radio data) and 2.8, and a break around a particle energy of 360 GeV, as required by the H.E.S.S. 
observations. 
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For a lower magnetic field in Zone 2 (20 — 30 /xG), the level of the observed 7-ray flux can 
be attributed mainly to this Zone. We cannot rule out either scenario with our current knowledge 
of the magnetic field. The shape of the TeV spectrum in our model could be improved with a 
combination of such relatively low magnetic field values in both Zones. 

The more interesting result obtained from our analysis of the observations is related to the 
constraints imposed by cumulative observations from the Fermi LAT, which allow us to con- 
clude that the TeV emission from SN 1006 cannot be produced by n° decay from a hadronic 
population with a "soft" power-law distribution (i.e., with an index ~ 2.0 or higher), such as 
the ones typically u sed in similar phenomenological mod e ls of broad-band SEDs of SNRs (e.g. , 
Arava & Cui 201o[ kbdo et al. 2010al : kbdo et al. 2010cl: kbdo et al. 2010el : krava etal. 2011 : 
Giordano et al. 2012|) . but they are consistent with a mixed leptonic-hadronic (or leptonic-dominated) 
model for the case where the index of the high-energy ion component is around 2, or with a pure 
leptonic model for softer cosmic ray distributions, as is seen in Figs. [2]and[3l Of course, the steeper 
the hadronic distribution, the lower its contribution to the TeV 7-ray flux should be. We cannot 
rule out a hard hadronic radiation spectrum (i.e., with a particle index below 2). 

It is important to stress the fact that our model only rules out hadronic distributions that are 
single steep power-laws at all energies as the main cause of the TeV emission from the source, and 
which differ from the ones predicted by some nonlinear models with efficient acceleration (e.g., 
Chevalier 1983c iMalkov 19991 : iBlondin & Ellison 200 1|) which show spectral hardening at high 
energies. H owever, such proton distributions are not seen in other high-efficiency acceleration 
models (e.g.. [Ellison et al. 2010|) . 



We also point out that our model is an extremely simple approximation which neglects all 
aspects of the temporal dynamics as well as nonlinear effec ts from the backreaction of accelerated 
particles in the fluid dynamics (see, e.g.. I Allen et al. 20081 for a description of spectral curvature 
seen in the radio spectrum of SN 1006 and possibly resulting from nonlinear dynamics), MHD 
turbulence and source photons other than the CMB that could be up-scattered and contribute to the 
IC flux. In fact, as can be seen in Fig. [TJ the shape of the H.E.S.S. spectrum of SN 1006 is hard to 
reproduce with our model, which is true for any other phenomenological IC-CMB emission model 
unless a rather odd spectral shape is used for electrons in Zone 1 (e.g., a distribution with a spectral 
break at low energies and a very steep spectrum, with a power- law index of ~ 3.7, above the break 
energy). A full broadband nonlinear model of the non-thermal emission could perhaps account for 
the spectral shape at TeV energies, but our main conclusion would still be valid. 

Based on our model it can be pointed out, as explained before, that a lower magnetic field 
in the X-ray emitting region is supported by the H.E.S.S. observations if the nature of the TeV 
emission is leptonic. The bipolar morphology of the 7-ray emission seen by H.E.S.S. indicates 
that it originates in the polar caps of the remnant, and in their proposed leptonic scenario, the flux 
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level is consistent with a field value of 30 /iG as pointed out by the authors (|Acero et al. 20101) . For 
their hadronic scenario, they find that a field of 120 /iG is consistent with the data, and necessary 
to "suppress" the leptonic emission. Our results from Section [3T21 are consistent with their leptonic 
emission scenario. We have already pointed out that a field value of 30 /iG in Zone 2 of our model 
can also account for the H.E.S.S. observations. 

On the other hand, other authors have derived a higher value for the magnetic field. Using a 
detailed nonlinear kinetic model for its non-ther mal X-ray emissio n, some estimates of the mag- 



netic field in SN 1006 favor a value of 150 /iG (IVolk et al. 20051) . which then translates into a 
predicted IC flux that is lower than the 7— ray flux observed by H.E.S.S., as can also be concluded 
with our model. This has lead ot her authors to argue that the TeV emission might be predominantly 



hadronic (|Berezhko et al. 2009|) . We can only point out that if such were the case, then our GeV 
limit on the emission implies that the hadronic 7-ray spectrum should be hard at GeV energies, 
which is our main co nclusion, and this is in accordance with the model shown by these authors 



( Berezhko et al. 2 009). 



The main leptonic scenario for the TeV emission in our two-Zone model presents other diffi- 
culties. The break in the particle spectra in Zone 1 which is necessary to reproduce the observations 
cannot be accounted for by synchrotron cooling alone. The particle break was introduced in the ba- 
sis of observations. A break in the predicted IC emission would be in general required by the LAT 
constraints and the break energy itself would be needed to reproduce the observed TeV spectral 
shape. 

With respect to the predicted bremsstrahlung flux, the fact that this emission (which repro- 
duces the particle spectrum directly) from SN 1006 is found to be negligible for all the leptonic 
parameters shown here means that the spectrum of the source at GeV energies would either be hard 
and leptonic in nature, hard and hadronic in nature or a relatively flat spectrum of mixed origin. The 
results of the Fermi LAT data analysis shown in Tabled] seem to indicate that the data favor a hard 
spectrum for the emission (as shown by the increasing values of the TS for lower photon indices) , 



which would also be consistent with the hadronic model mentioned above (|Berezhko et al. 20 09). 
This will only be confirmed with a significant GeV detection of the source in the future. 

With respect to hadronic emission from a particle population with an index of 2 (Fig. [3), it 
is seen that the data are more consistent with mixed leptonic and hadronic contributions (e.g., of 
comparable fluxes at TeV energies) rather than a predominantly hadronic scenario, for which the 
resulting fluxes are slightly above the GeV upper limit. 

Finally, we point out that despite its problems, a leptonic origin of the non-thermal emission 
from SN 1006 at all ene rgies would be consis tent with the fact that other "SN 1006-like" SNRs 



such as RX J17 13.7-3946 (|Koyama et al. 19971) . show 7-ray spectra and shock structure consistent 
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with leptonic emission (|Ellison et al. 2010uAbdo et al. 201 la|) . Also, it has been proposed recently 
that non-thermal emission from S NRs in low dens ity environments (such as SN 1006) is caused 
primarily by high-energy leptons (|Yuan et al. 2012|) . 



We cannot of course conclude that the TeV emission from the source is mainly hadronic or 
leptonic, but we stress the fact that the data rule out certain hadronic emission models such as the 
ones that have been used to account f or the observations of other young SNR s, such as Cas A 
jAbdo et al. 2010bl : lArava & Cui 201Clh and Tycho's SNR dGiordano et al. 20121) . 
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Table 1: Fermi LAT Upper Limits on the GeV SED of SN 1006 for Different Assumed Photon 
Spectral Indices 



Index 


£^ UL (ergs cm- 2 s" 1 ) 


TS 


1.6 


1.0 x 10~ 12 a 


7.2 


2.0 


9.6 x 1(T 13 


5.7 


2.14 


9.9 x 10~ 13 b 


4.9 



a Evaluated at an energy of 5.7 GeV. 
b Evaluated at an energy of 2. 1 GeV. 



-17- 



10"'° T 




Fig. 1. — SED of SN 1006 and leptonic emission model. The model includes synchrotron and 
IC/CMB emission from Zone 1 (dashed lines) and from Zone 2 (dotted lines), and total emission 
(solid line). The parameters of the fit are reported in Sections [3.11 and [3T2T . 
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Fig. 2. — SED of SN 1006 and models with a contribution from a hadronic component produced by 
a power-law distribution of particles with index 2.14 (dash-dotted line), (a) assuming a hadronic 
origin for the TeV emission; and (b) showing the maximum hadronic 7-ray flux allowed by the 
data. 
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Fig. 3. — Same as Fig. [2] for models with a hadronic component produced by an ion distribution 
with index 2.0 (dash-dotted line), corresponding to (a) the hadronic-dominated scenario; and (b) 
the mixed leptonic-hadronic case. 



